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ABSTRACT 

Context. The 3-torus (r 3 ) Friedmann-Lemaitre-Robertson-Walker (FLRW) model better fits the nearly zero large-scale two-point 
auto-correlation of the Wilkinson Microwave Anisotropy Probe (WMAP) cosmic microwave background (CMB) sky maps than the 
infinite flat model. The r 3 model's parameters, recently found using an optimal cross-correlation method on WMAP data, imply 
approximately equal-redshift topological lensing at redshifts z ~ 6, the redshift range of the upcoming generation of new instruments 
and telescopes. 

Aims. We investigate observational strategies that can reject the T 3 solution for a given region of parameter space of physical 
assumptions, or provide good candidate topologically lensed galaxy pairs for detailed spectroscopic followup. 

Methods. T 3 holonomies are applied to (i) existing z ~ 6 observations and (ii) simulated observations, creating multiply connected 
catalogues. Corresponding simply connected catalogues are generated. The simulated observational strategies are motivated by the 
matched discs principle. Each catalogue is analysed using a successive filter method and collecting matched quadruples. Quadruple 
statistics between the multiply and simply connected catalogues are compared. 

Results. The expected rejection of the hypothesis, or detection of candidate topologically lensed galaxies, is possible at a significance 
of 5% for a pair of T 3 axis-centred northern and southern surveys if photometric redshift accuracy is cr(zphot) ~ 0.01 for a pair of nearly 
complete 100 deg 2 surveys with a total of * 500 galaxies over 4.3 < z < 6.6, or for a pair of 196 deg 2 surveys with > 400 galaxies 
and 0"(Zphot) ~ 0.02 over 4 < z < 7. Dropping the maximum time interval in a pair from At = 1 h Gyr to At = 0. 1 h~ l Gyr requires 
cfephot) ~ 0.005 or c(z phot ) < 0.01, respectively. 

Conclusions. Millions of z ~ 6 galaxies will be observed over fields of these sizes during the coming decades, implying much 
stronger constraints. The question is not ;/the hypothesis will be rejected or confirmed, it is when. 

Key words, cosmology: observations - cosmological parameters - Galaxies: high-redshift — distance scale 
1. Introduction Spergel et al. I2003I also Fig. 1 Roukema et al. I2008bl i where 
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■ rsLS is the comoving distance to the surface of last scattering 
The predictions implied by cosmic topology interpretations of (SLS) xhuSi Occam's razor suggests that an approximately zero 
T the lack of large-scale (hereafter, r > lO/^ 1 Gpc) structure in large . scale auto . colTe lation function provides a simpler statisti- 
cs) the cosmic microwave background (CMB) maps are relevant cal model for the largest scalefj than the harmonic ana i ysis . 
© for the design of deep redshift observational strategies. The The Friedmann-Lemaitre-Robertson- Walker (FLRW) mod- 
Cn large-scale structures in the CMB as observed by the COsmic elfj of ±e Universe assume that the comoving spatial sect ion 
Background Explorer (COBE) and the Wilkmson Microwave ifj a 3 _ manifold> having both curvature an d topology (e.g., de 
> Amsoti-opy Probe (WMAP) (Bennett et al. [2003}, i.e., in par- shter rjgjTJ priedmann [J923] [USS Lemaitre [19271 Robert- 
'O tlcular ' the second moment of the tem P erature fluctuation distn- SQn rr^f jTl For recent reyiews and terminology such as « fun _ 
K bution, can be statistically analysed using either (1) spherical har- damental domain ^ << cove ring space", and "topological lensing" 
H monies analysis of the temperature fluctuations, which is invalid (where the topology of a spatial section of the UniverS e can be 
if observations near the galactic plane are contaminated, or (11) thought of afj a fa analogy wkh gravitational lensing)? see 
the angular or spatial two-point auto-correlation function of the Lachieze-Rey & Luminet (TT995T> : Luminet (TT99ST> : Stark- 
temperature fluctuations. (The full moment hierarchy can be es- man (Tjgggj. Lumine t & Roukema ( [19991 : Rebou e as & Gomero 
timated using Minkowski functionalsof iso-temperature excur- Qr shorter> Roukema (goooj. Luminet mm . 
sion sets (Schmalzing & Buchert [19971 Schmalzing & Gorski structures bigger than the fundamental domain of an FLRW 
mm Ducout et al. [Ml.) Interpreting the (second moment) model (i e of the 3 _ manifold of the comoving spatia i sec . 
large-scale harmonics to be generated by statistically mdepen- tion) cannQt physically exist) so observa tions that are most eas- 
dent Gaussian distributions requires an unlikely '^piracy" Uy interpreted m the universa i covering space (apparent space) 
between the harmonics at low I values(Copi et al. [20071 [2009J should approximate i y revea l this (Starobinsky H993] Stevens 
Sarkar et al. [2011] Copi et al. [2010) in order to match the et d rj^ though see also Fig h Roukema [J996>. Thus, 
nearly zero value of the (two-point) auto-correlation function at 

large (r > r SLS * 10.0/i _1 Gpc) pair separations (Sect. 7, Fig. 16 1 See Lemaitre dT93Tl i for an English translation of Lemaltre (fT927l . 

with Lemaitre's observational estimate H » 600km/s/Mpc excluded 

* BFR: during visiting lectureship. (e.g. van den Bergh 120111 Block [2"0TTI l. 
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Table 1. Preferred fundamental directions 1 < i < 6 of the Aurich 
d2008t T 3 solution in galactic coordinates 3 , with a directional uncer- 
tainty of ~ 2° (great circle degrees). 
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unit 
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3 4* 


5* 


6 h 
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17 


107 186 


197 


287 
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77 


-13 


3 -77 


13 


-3 


Z/ 3 


c/H 






3.80 + 0.05 c 







Notes. w The typographical error (sign of b T 6 ) in Table 1 of Aurich 

d2008t has been corrected (cf Figs 6, 7 of Aurich 120081 

<b) The directions 4, 5, and 6 are exactly antipodal to the directions 1, 2, 

and 3, respectively, and are listed for convenience only. 

(c) Mean and standard error adopted from our analysis (12.1 1 . 

the simpler interpretation of the large-scale CMB observations 
is that the Universe is a multiply connected FLRW model rather 
than a simply connected FLRW model: the near-zero large- 
scale two-point auto-correlation reveals the size of the Universe. 
Among the models that are considered to provide better fits to the 
WMAP data than the infinite flat model are the Poincare dodec- 
ahedra l space {S 3 /I* Lum inet et al. 12003] Auri ch et al. 12005 al 
I2005bl Gunde rmann 120031 Caillerie et al. 120071 Roukema etal. 
I2008bl l2008aT>, the Picard space (// 3 /PSL 2 (Z[/]), Picard [T884l 
Aurich et al. |20"04| |, the re gular 3 -torus (T 3 , S pergel et al. 120031 
Aurich et al. 120071 Aurich |2008l Aurich et al. EoTjg] EoTUt . and 
the "2-torus" (S'xS'xR, Aslanyan & Manohar l2012l i, though 
some analyses, also assuming an FLRW (homogeneous) metric, 
favour the infinite flat model (Key et al. 120071 Niarchou & Jaffe 
120071 Bielewicz & Banday |20TTT i. 

While there is a lack of consensus in fitting models to the 
WMAP data, the T 3 model fits have a particular advantage from 
an observational point of view: they imply sub-SLS topologi- 
cal lensing at redshifts that are becoming observationally realis- 
tic. Many cosmic topology CMB analyses consider ensembles 
of possible sizes and orientations of the fundamental domain of 
a given model. A more empirical approach requires that the fun- 
damental domain of the real Universe must have a specific size 
and orientation in astronomical coordinates. We live in a specific 
realisation of the physical processes that led to our Universe, 
not an ensemble of realisations. Aurich (12008) applied the opti- 
mal cross-correlation method (Roukema et al. I2008bl 12008a) to 
find a specific T 3 solution. This solution implies equal-redshift 
topological lensing at redshifts z i 6, the redshift range at which 
many detections are expected with the upcoming generation of 
new instruments and telescopes. Thus, this candidate 3-manifold 
is, in principle, empirically testable independently of the CMB. 

Roukema & Kazimierczak ( 1201 II ) presented a corollary of 
the matched circles principle (Cornish et al. 119961 11998) : the 
matched discs principle. In general, topologically lensed pairs 
of images of a given physical object occur at different redshifts. 
This complicates observational tests: the lifetimes of quasars are 
short compared to typical differences in lookback times, and at 
high redshifts, a galaxy seen in one direction at a redshift of 
Z\ may be absent at an expected position with a higher red- 
shift zi > Z\ because the initial starburst occurs after zi and 
before z\- Matched discs minimise this problem, by selecting 
the set of spatial positions — discs — for which a topologically 
lensed pair of any given galaxy occurs at identical redshifts in 
the two matched discs (Fig. 1, Roukema & Kazimierczak l20 1 Tt . 



The redshift is lowest at the centre of a disc and increases radi- 
ally outwards. The Poincare dodecahedral space matched discs 
occur with z m \ n = 106 + 18, so detection of gravitationally col- 
lapsed objects (galaxies) would require the existence and detec- 
tion of rare high overdensity peaks of the density fluctuation 
distribution that collapsed very early. In contrast, the Aurich 
( 120081 T 3 candidate has z min ~ 5 (Fig. 10, Sect. 3, Aurich[2008]l. 
While relatively small numbers of quasars have been detected 
with z i 5, many Lyman break galaxies (LBGs) and Lyman alpha 
emitters (LAEs) in this redshift range have now been observed 
[with spectroscopic confirmations up to z = 8.6 (Lehnert et al. 
120 10H and many more are likely to be detected with the new in- 
struments and telescopes of the coming decade. LBGs and LAEs 
have a topological lensing advantage over quasars: the emission 
from a given object is expected to be much more isotropic than 
the beamed jets of a quasar, and the lifetime is likely to be much 
longer, because of long galaxy dynamical time scales and stellar 
lifetimes. 

Thus, the Aurich ( 12008b solution is going to be testable over 
matched discs of many square degrees using observations to 
be made over the next few years (e.g., VISTA/VIKING, Sub- 
aru/ CISCO , EUCLID, VLT/X-Shooter, VLT/MUSE, JWST; see 
Sect. l2.2.2l for details). Although surveys that are approximately 
limited to a narrow band in redshift correspond to thin shells 
rather than thin discs, surveys over solid angles <k 4tt (e.g. 
many square degrees) can strongly overlap with matched discs 
if made over the predicted redshift ranges. Since observational 
and theoretical uncertainties in the candidate 3-manifold cannot 
be avoided, a successive filter method for finding sets of likely 
topologically lensed pairs (Roukema 119961 Uzan et al. 119991 
Marecki et al. 120051 Fujii & Yoshii 1201 Ti l is preferabl e to the 
pair separation histogram (PSH) method (Lehoucq et al. 11996) , 
as the former can reveal very weak topological lensing signals 
(Fujii & Yoshii 1201 1) . Here, a successive filter method, moti- 
vated by the matched discs principle, is applied to existing ob- 
servational data and planned or possible surveys, using numer- 
ical simulations, in order to investigate observational strategies 
that can reject the T 3 solution, on the assumption that the FLRW 
metric is close enough to physical reality. 

In Sect. |2J] we briefly present the Aurich (12008) T 3 so- 
lution, matched discs, and not-quite-matched beams. Existing 
and planned or possible observations that can be used to pre- 
dict sky positions and redshifts of topologically lensed copies 
are presented in Sect. 12.21 For existing observations, we apply 
the successive filter method to two partly simulated catalogues. 
The signal is represented by analysing the union of the observa- 
tional catalogue with its topologically lensed images on (approx- 
imately) the opposite side(s) of the sky. Noise is represented by 
analysing the union of the observational catalogue with a sim- 
ulated (non-lensed) data set in the region of comoving space 
where lensed objects would be found. Thus, signal is compared 
to noise. For planned or possible observations, the data are fully 
simulated, comparing a multiply connected T 3 simulation with a 
simply connected R 3 simulation. The simulation methods are de- 
scribed in Sect. 12.31 and the successive filter method is presented 
in Sect. 12.41 Statistical tests, either to find a low probability of 
falsely excluding the T 3 hypothesis, or to find a low probabil- 
ity of falsely detecting evidence in favour of the hypothesis, are 
presented in Sect. 12.51 

Results are presented in Sect. |3]and conclusions are given in 
Sect. |4] All distances are FLRW comoving distances except if 
stated otherwise, and Q m and Qa are the dimensionless matter 
density and dark energy density, respectively. The Hubble con- 
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Fig. 1. Sub-gigaparsec-scale mean cross-correlation (fc>,- (L) 

[Eq. (T}] in pK 2 as a function of Universe size L T ' in units of c/Hq, 
for the WMAP 9-year ILC map using the KQ85 galactic mask, for 
r = 0.2,0.4, 1.0 h~ l Gpc from top to bottom (red, green, blue, respec- 
tively, online), for the orientation given in Table[T] 




matched discs 



Fig. 2. Matched discs in the universal covering space, as per Fig. 1, 
Roukema & Kazimierczak (2011b . but for the Aurich d2008t T 3 solu- 
tion, with the addition of matched pencil beam observations. Spheres 
corresponding to an example redshift z = 6.2 are shown intersecting 
themselves (in the covering space). The redshift at the centres of the 
matched discs is less than this, and within a disc, the redshift increases 
radially out to the surface of last scattering (SLS). An observer (at the 
centre) pointing a telescope towards the right-hand matched disc can 
observe the z i 6.2 portion of the right-hand copy of the pencil beam at 
roughly the same redshifts as those of the corresponding portion of the 
left-hand copy of the pencil beam. 



stant is written Hq = lOOh km/s/Mpc and c is the conversion 
factor from time units to space units. 



Table 2. Redshifts increasing radially outwards from the centre of 
matched discs for the Aurich J2008) T 3 solution, but with L T = 3.8 
motivated by the WMAP7 ILC map (Fig. [T}, at angle j} from the disc 
centre, and approximate fraction of full sky covered by the three pairs 
of discs LL>/(4n) a , depending on the assumed FLRW metric parameter 
fi m , where £1 A = 1 - £2 m . 



O m — 


0.26 


0.27 


0.28 


0.29 


0.30 




p 






z 






w/(4tt) 


o 





4.98 


5.17 


5.35 


5.55 


5.75 


0% 


5 


5.04 


5.22 


5.42 


5.61 


5.82 


1% 


10 


5.21 


5.41 


5.61 


5.82 


6.04 


5% 


1& 


6.10 


6.36 


6.63 


6.92 


7.27 


5% 


15 


5.52 


5.74 


5.96 


6.20 


6.45 


10% 


20 


6.02 


6.27 


6.53 


6.81 


7.10 


18% 


25 


6.78 


7.1 


7.42 


7.77 


8.13 


28% 


30 


7.99 


8.41 


8.85 


9.31 


9.81 


40% 



Notes. (a) Galactic masking would modify these rough estimates. 
<b) Calculation assuming L T = 4.0 for the circle with p = 10°, to illus- 
trate the effect of uncertainty in L T . 

2. Method 

2.1. The Aurich \2008\ T 3 solution, matched discs, and 
not-quite-matched beams 

The coordinates of the Aurich (I2008l > T 3 solution are given in Ta- 
bleQ] We check the cross-correlation £c (Roukema et al. I2008b| ) 
on the 9-year WMAP internal linear combination (ILC) mapQ, 
using the corresponding KQ85 galactic contamination mask, 
and show the dependence of the sub-gigaparsec-scale mean 
cross-correlation 

m r {L):= - f tc{r',L)&r' (1) 
r Jo 

in Fig. Q] The cross-correlation £c is a two-point correlation 
function that at small scales correlates pairs of points that are 
observed to be separated by large distances but according to 
the hypothesised 3-manifold are also separated by a small dis- 
tance. The cross-correlation should be low if the hypothesis is 
wrong (Fig. 2 Roukema et al. I2008bb and high if the hypoth- 
esis is correct (Fig. 3 Roukema et al. I2008bb . The strongest 
cross-correlation in Fig. [T]occurs for L T « 3.80 + 0.05 (in units 
of c/Ho), where the uncertainty is the half-width of the maxi- 
mum at its base. The Aurich ( 12008b uncertainties are estimated 
from the pseudo-probability function used in the Monte Carlo 
Markov Chain search for an optimal solution. Since the pseudo- 
probability estimator is not a true probability, other methods are 
needed to estimate a realistic uncertainty. Figure 9 of Aurich 
( 12008b . showing the estimates in different foreground-subtracted 
single waveband WMAP maps, suggests cr(L ) ~ 0.02, com- 
bining random and some systematic sources of error (ILC ver- 
sus single band maps). The shift by 0.05(c/Hq) between Au- 
rich ( I2008b 's estimate of L 7 " 3 for WMAP5 data and that shown 

2 http : / /lambda . gsf c . nasa . gov/data/map/ dr 5/dfp/ 
ilc/wmap_ilc_9yr_v5 . fits 

3 The 9-year version of the KQ85 mask 
http://lambda.gsfc.nasa.gov/data/map/dr5/ancillary/masks/ 
wmap_temperature_kq85_analysis_mask_r9_9yr_v5 . fits 

hides about 25.2% of the sky. 
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in Fig. Q] suggests that inclusion of systematic error would give 
cr(L r3 ) > 0.05, i.e. 150/T 1 Mpc, rath er than o-(L r3 ) = 0.02. The 
angular uncertainty given by Aurich (12008) . 2°, is similar to that 
of the Poincare dodecahedral space solution found earlier with 
the optimal cross-correlation method (Roukema et al. I2008bl 
I2008al ). i.e. a few hundred comoving hr x Mpc in a tangential 
direction at a radial distance of 5 to 10 h~ l Gpc. 

Leaving aside these uncertainties for the moment, the obser- 
vationally most dramatic element of this T 3 solution is revealed 
by considering the redshift of an object topologically lensed in 
opposite directions along a fundamental axis. Aurich (12008b 
shows this in his Fig. 10 and briefly discusses it. Here, we use 
the principle of matched discs (Roukema & Kazimierczak l201 il l, 
shown with some extra detail in Fig. [2] and Table [2] listing the 
redshifts at the centres of a pair of matched discs and at circles 
increasing radially outwards to form the pair of matched discs. 
The fractional sky coverage <y/(47r) = 12^(1 - coSj6)/(47r), for 
discs of angular radius ft, is also shown. 

Given the rapidly increasing numbers of observed objects in 
the 5 < z ~ 6 range, the beginning of detections around 6 < z 8, 
and the large fractional covering of the sky by the matched discs, 
it is obvious that independently of the few hundred hr x Mpc 
uncertainties in the T 3 solution, the observational data to re- 
ject or confirm the solution to very high statistical significance 
are going to become available over the next decade or two. 
The absence of topologically lensed images of quasars can be 
interpreted as a problem of their short lifetimes and highly 
anisotropic nature (beamed jets), but the absence of topologi- 
cally lensed images of early forming galaxies would be much 
more difficult to explain. 

Let us return to the uncertainty in the T 3 solution parameters. 
Pencil beam observations ("deep fields") that probe to z ~ 6 have 
typical widths of at most a few arcminutes, i.e. a few hr x Mpc 
in comoving thickness at these redshifts. Thus, an accuracy of a 
few hundred h^ 1 Mpc suggests that deep fields of a few degrees 
in size — pointed in the appropriate directions — will most likely 
be needed. The fourth row of entries in Table |2] shows the red- 
shifts at the f3 = 10° radius circle for L — 4.0, i.e. 4<x greater 
than the WMAP9-ILC estimate shown in Fig. [T] if cr is defined 
as the absolute difference between the Aurich (12008) WMAP5 
estimate and the WMAP9-ILC estimate made here. Thus, for 
surveys limited to, e.g. z < 6.2, observations closer to the face 
centres would better test the solution. 

Figure [2] also shows a pair of not-quite-matched beams. 
These are not truly matched, since they exactly match only 
where they intersect the matched discs. The topologically 
matched parts of the two beams to the left (in Fig. |2]i of their re- 
spective discs occur in the left-hand copy of the beam at slightly 
higher redshifts and in the right-hand copy at slightly lower red- 
shifts, with respect to the redshift of the intersections with the 
matched discs. Since realistic surveys usually have wide redshift 
distributions, there is a fair chance of covering the topologically 
lensed region. However, this requires that the survey of the right- 
hand copy (in this figure) of the beam be wide enough in solid 
angle to cover the projection of the "beam" onto the sky, since 
the right-hand copy is not a beam from the observer's point of 
view. For a small angular distance of the left-hand copy of the 
beam from the matched disc centre, the projected solid angle 
of a small portion of the right-hand copy will not be too large. 
Thus, for a small angular offset from the matched disc centres 
(axes of the fundamental domain), "not-quite-matched beams" 
can potentially test topological lensing hypotheses. 

While matched discs and not-quite-matched beams indicate 
the parts of the sky that should be observed, detecting a signif- 



Table 3. Number of astrophysical objects with 5 < z < 6.2 and within 
10° of the Aurich ( 12008) T 3 solution fundamental axes (numbered as in 
Tableffl listed in NED" as of 28 Oct 2012. 



i 1 2 


3 4 


5 6 


NED count 44 


9 


1 



(a) 



Notes. 

( http : / /ned .ipac.caltech. edu ) 



NASA/IPAC Extragalactic 



Database, 



Table 4. Some well-known deep fields and their angular distance 
from the nearest T 3 axis (listed if < 20°). 



name 


I 


b 


min^Ce,) 


If 3 




unit 


o 


o 


o 


o 


o 


VVDS 0226-04" 


172. 


-58. 


20. 


186. 


-11. 


CDFS* 


224. 


-54. 


26. 






HDF-S C 


328. 


-49. 


52. 






VVDS 2217+00" 


63. 


-44. 


50. 






VVDS 1003+01" 


238. 


43. 


46. 






HDF J 


126. 


55. 


43. 






AEGIS f 


97. 


60. 


33. 






VVDS 1400+05 / 


343. 


63. 


16. 


6. 


11. 


SDF ff 


38. 


83. 


8. 


6. 


11. 



Notes. {a) VIRMOS Very Deep Survey (Le Fevre et al. 120031 

(b) Chandra Deep Field South 

(http: //www. eso.org/~vmainier/cdfs_pub/CDFSfield.html ) 

(c) Hubble Deep Field-South "WPFC2 (apex) pointing" 
(http : //www. stsci .edu/ftp/science/hdfsouth/coordinatesS.html I 
<d) Hubble Deep Field (Williams et al. [19961 

<e) All- Wavelength Extended Groth Strip International Survey (Davis 

etal. 120071 1 

<f) Iovino et al. ( l2005l > 

(g) Subaru Deep Field (Maihara et al. l200Tt 



icant statistical signal and generating a list of candidate pairs 
of topologically lensed objects requires a sensitive statistical 
method. This is presented below in Sect. 12.41 However, first 
we need to consider existing, planned and possible surveys in 
the three-dimensional regions of interest. 



2.2. Observational catalogues 
2.2.1. Existing observations 

Table [3] shows that almost all of the z ~ 6 objects that are so 
far known near the T 3 axes are those near the high galactic lati- 
tude, northern fundamental direction 1 , although a few are also 
known near the corresponding southern fundamental direction 
4. Most of the northern objects are from the Subaru Deep Field 
observed with the Cooled Infrared Spectrograph and Camera 
for [the] OH airglow suppressor (CISCO) on the 8.2 m Subaru 
telescope (Maihara et al. 120011 Shimasaku et al. 12006b . listed 
along with some of the other well-known deep fields in Table [4] 
The galactic coordinates and redshifts of the northern objects are 
listed as part of Table [5] 
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Table 5. Known objects" near the northern galactic matched disc at sky position (h,bi) and redshift zi and their predicted most likely southern 
galactic positions (fe, b^) and redshifts zi, and redshift and cosmological time differences Az := zi - Zi and Af := t2 - t\ (in h~ x Myr), respectively, 
assuming metric parameters fi m = 0.28, Q r = 1.65 x 10~ 4 fl m ,fl A = 1 - f2 m - fl r and the T 3 fundamental axis (/ = 6°,b = +77°) of length 
L = 11.4/r' Gpc. 



name 


h 


bi 


h 


bi 




Z\ 




Z2 




A - 


At 


[SKD2006] 017721 


35.597 


82 


.506 


174 


.850 


-70 


.004 


5 


.724 


5. 


,281 





,443 


75 


.6 


[SKD2006] 036334 


36.053 


82 


.455 


174 


.632 


-70 


.020 


5 


.705 


5. 


,300 





.405 


69 


.2 


[SKD2006] 034503 


36.000 


82 


.515 


174 


.771 


-69 


.984 


5 


.698 


5 


.308 





.390 


66 


.6 


SDFJ132519.4+271829 


36.238 


82 


.470 


174 


.368 


-69 


.758 


6 


.000 


5 


.058 





.942 


160 


.2 


[SKD2006] 018699 


35.605 


82 


.619 


175 


.130 


-69 


.980 


5 


.645 


5 


.356 





.289 


49 


.4 


[SKD2006] 020087 


35.642 


82 


.613 


175 


.045 


-69 


.917 


5 


.724 


5. 


,286 





,438 


74 


.8 


[SKD2006] 031858 


35.934 


82 


.554 


174 


.870 


-69 


.965 


5 


.692 


5. 


,314 





,378 


64 


.5 


SDSS J 1 34040.24+28 1 328.1* 


41.523 


79 


.054 


165 


.545 


-71 


.950 


5 


.338 


5. 


,654 


-0 


.315 


-53 


.9 


[SKD2006] 020495 


35.635 


82 


.719 


175 


.272 


-69 


.858 


5 


.697 


5. 


,314 





.383 


65 


.4 


SDFJ132426.5+271600 


35.900 


82 


.666 
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Notes. (a> The objects are galaxies except where otherwise indicated. <b) QSO. (c) Object type uncertain. 



2.2.2. Planned and possible observations 

First let us consider high-redshift quasars. The VISTA/VIKING 
4-m class telescope project should discover quasars at z ~ 7 over 
about 1500 deg 2 centred on the northern and southern Galac- 
tic polesQ (Findlay et al. 120111 1. Although Table |2] here and 



4 The criterion "centred on the northern Galactic cap" presumably 
means "a large solid angle centred on the northern Galactic pole". 



the expected completeness levels shown in Fig. 13 of Findlay 
et al. (1201 II ) indicate that coverage of the desired regions is pos- 
sible (depending on the value of Q. m ), only about 8 quasars are 
expected in the whole survey (Sect. 4.1, Findlay et al. 1201 U . 
clearly too low to obtain significant evidence either for or against 
the T candidate. 

Polsterer et al. (120121 1 find 22,992 photometric quasar can- 
didates with 5.5 < z < 6.2 and an estimated error of cr(Az) = 
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0.087, expecting about half the candidates to be true quasars. 
However, the solid angle from which the quasars are selected is 
about 2tt, giving a candidate surface density of about 1/deg 2 and 
an expected spectroscopic number density of about 0.5/deg 2 if 
a complete spectroscopic followup were to be performed. This 
is about 100 times higher than for the VISTA/VIKING survey 
(as presently designed), but still somewhat low. While quasars 
dominated high-redshift records for the decades when z > 1 was 
considered a high redshift for a spectroscopically confirmed ex- 
tragalactic object, this seems less likely at z ~ 6. For a given ap- 
parent magnitude limit, it is clear that LBG and LAE catalogues 
are present in much high number densities at z ~ 6. 

Plans for LBG and LAE searches typically focus on the ex- 
isting "deep fields", with the aims of achieving wide coverage 
across the electromagnetic spectrum. Several of these fields and 
their angular distances from the closest T 3 axes are given in 
Table [4] indicating which fields would be the most useful for 
testing the T 3 candidate via LBG and LAE searches. Compar- 
ison with Table [2] indicates that for 5 < z i 7 searches, the Sub- 
aru Deep Field and its corresponding northern region is optimal, 
while for 5.5 i 8 searches, the Virmos Very Deep Survey fields 
0226-04 and 1400+05 and their northern and southern, respec- 
tively, counterparts would be worth observing. 

The EUCLID mission^] planned for launch in 2O1S0 has a 
deep survey (Sect. 1.3, Refregier et al. I20T01 the EUCLID 
optical and NIR Deep Imaging Survey, planned over 40 deg 2 . 
This should obtain "thousands" of likely LBG's and LAE's with 
z > 6 based on photometric redshifts (Sect. 14.2, Refregier et al. 
120101 1. i.e., >50/deg 2 . If a few objects had quadruple multi- 
plicities that made them very likely topological lensing candi- 
dates, then optical/near-infrared spectroscopic followup with an 
instrument such as X-Shooter (300 nm < A < 2400 nm, Vemet 
et al. 1201 lb on the Very Large Telescope (VLT) on the southern 
objects and corresponding spectra of the northern objects with 
the Keck or Subaru telescopes should obtain enough rest-frame 
spectral energy distribution information to check whether the 
would-be topologically lensed pairs of objects resembled each 
other more than would be expected for arbitrary pairs of ob- 
jects at similar redshifts. The Multi Unit Spectroscopic Explorer 
(MUSEd on the VLT, which should be well-tested by the time 
that the EUCLID data are available, would be useful for sec- 
ondary followup to study the spatial environments of the candi- 
date lensed objects, especially when the cosmological time dif- 
ference between the members of a pair is small compared to a 
typical galaxy dynamical time scale. 

Thus, we should simulate several square degrees of obser- 
vations with parameters consistent with those estimated for EU- 
CLID. In particular, it would be interesting to see if a statistical 
signal could be obtained from photometric redshifts alone, prior 
to spending high amounts of exposure time on highly sought af- 
ter telescope/instrument combinations. 

The Ultra-Deep Survey of the James Webb Space Telescope 
(JWST) should detect galaxies at z i 6, but over only about 
10 ai-cmin 2 (Sect. 2, Table II, Gardner et al. l2006l >. The JWST's 
planned Deep Wide Survey should find galaxies over 1 < z i 6 
over a larger solid angle, 100 arcmin 2 (Sect. 3.6, Table III, Gard- 
ner et al. 120061 1. For an initial uncertainty of about two degrees 
in the T 3 axes, these surveys are clearly too narrow. 



5 http://sci.esa.int/eucli5 

6 http://sci . esa. int/science-e/www/ object/ 



2.3. Simulations 

Simulating searches for type I quadruples requires comparison 
of observations in a T 3 model to those in an R 3 model. The ex- 
isting Subaru Deep Field observations are used to provide an ob- 
servationally based simulation of both sorts, in which the south- 
ern field is generated for the T 3 model by applying the Aurich 
(12008b T 3 holonomy (Table [TJ to the Subaru galaxies, and sim- 
ulated independently of real observations in the case of the K 3 
model. For a hypothetical observational program aimed at the 
centres of both the 1-4 axis expected matched discs, fully simu- 
lated data are needed in both cases. 

The spatial two-point auto-correlation function f (r) could, in 
principle, lead to excess non-topological chance isometries if the 
tolerances for requiring isometry are not tight enough. For the 
low number densities of objects expected, the effect is unlikely 
to be strong. Half of a large number of points are first generated 
uniformly within a redshift shell defined by the required redshift 
limits. Each further point is chosen by randomly choosing an 
existing point, and placing the new point in a 3-Euclidean di- 
rection chosen uniformly from An ster at a comoving distance 
chosen with the probability P(r, dr) oc 1 + g(r), where 



(?)'■ 

I 0, 



r < r„ 



<r<r„ 



(2) 



r > r„ 



for a correlation length ro and power law index y, and numerical 
cutoffs r mm and r max . Points that do not fall in the redshift shell 
are ignored, and new points are generated in the same way until 
the required number fill the shell. Numerical measurement of 
the resulting simulated distributions shows that in practice, this 
simulates a stronger correlation than that required. Thus, the 
chance of non-topological isometries is overestimated, i.e. our 
results are conservative: real observations are less likely to give 
a false positive detection. Moreover, in order to err on the side of 
possibly underestimating the numbers of matching quadruples in 
the multiply connected case, the correlation function is not used 
when simulating this case; uniform distributions are drawn from 
instead. Again, this is conservative: we slightly underestimate 
the statistical significance of the method (Sect. l23l l. 

Peculiar velocities and the uncertainties in redshift estima- 
tion need to be simulated (Roukema l 19961 ). Photometric redshift 
errors are those of most interest here: are they small enough to 
significantly discriminate between signal and noise? These er- 
rors are all simulated radially using a peculiar rapidity <p pec se- 
lected from a Gaussian distribution of mean zero and standard 
deviation atanh(cr(y8p ec )) for a single input parameter, the pe- 
culiar velocity standard deviation o-(J3 pec ), and using Eq. ( lA.4b . 
where y6p ec = tannic andz pec = [(l+ySp ec )/(l-y8pec)] 1/2 -l-Use 
of rapidities rather than velocities avoids the unphysical case of 
^8ped ^ 1 for simulating photometric redshifts with high <x(/? pec ). 

Before applying the holonomies in a given simulation, ran- 
dom Gaussian errors are added to the T 3 axis parameters given 
in Table Q] "Observations" of the simulated catalogues are car- 
ried out in galactic coordinate limited regions in the expected 
position(s). 



2.4. The successive filter method for type I pairs 



index. cfm?fobjectid=49385 
1 [Ettp : //muse . univ- lyonl . fr/ 



Lehoucq et al. ( 11999b introduced the terminology of type II pairs 
and type I pairs (or «-tuples), that had been introduced earlier by 
Lehoucq et al. ( 1 19961 ) and Roukema ( 11996b , respectively. Type 
I pairs or «-tuples can be thought of as the matching between a 
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local region of objects and its distant copy. Since a holonomy 
for an FLRW multiply connected model is an isometry, the dis- 
tances among the members within the "original" region should 
correspond to the distances among the members of the copy of 
that region. A Type II pair — in a space such as T 3 — corresponds 
to an object and its copy. In T 3 , this separation is a vector in 
the covering space R 3 . Difficulties in finding statistically sig- 
nificant numbers of either type of pair given realistic observa- 
tional parameters and noise led to a new way of collecting type I 
pairs (Uzan et al. 11999b and a successive filter method (Marecki 
et al. 120051 ). Fujii & Yoshii ( 1201 II ) presented a variation on 
the successive filter method, along with an analysis step that is 
roughly equivalent to collecting together n-tuples of mixed type 
for all n > 4 or to the Marecki et al. (120051 bunches of pairs 
filter (Sect. 3.2.5. Fuiii & Yoshii 1201 U . The latter step, made 
possible by the successive filters, significantly reduces the com- 
binatorial problem presented in Roukema (11996) . A catalogue 
containing two distant regions, each with N objects, might con- 
tain just one pair of matched n-tuples across the two regions, but 
searching for all ( N C„) 2 combinations of objects with N = 100 
and n — 7 would require 10 20 comparisons of 7-tuples. Fujii 
& Yoshii ( 1201 11 ) use toy simulations to show that a very small 
number of matched n-tuples can be detected by the full method 
of successive filtering and collecting. 

The successive filter method implemented here tests four 
simulated (and/or real) galaxies at (x, y, z)(i), (x,y,z)(j), 
(x, y, z)(k), and (x, y, z)(l) via the following filters in the follow- 
ing order, where (i, j, k, I) is an ordered choice of unequal indices 
in the list of galaxies. A crossed quadruple parameter ^f(q) e Z2 
is defined for a given ordered quadruple q = (z, j, k, I) in order 
to relate the lifetime and bunches of pairs (BoP) filters (Marecki 
et al. 120051 ). This algorithm allows the pair separations d, where 
d(., .) is the comoving distance between two arbitrary points in 
the covering space R 3 , the (signed) x, y,z component separa- 
tions of each pair, to be calculated first, and the loop for finding 
quadruples to be performed over pairs of pairs. Arithmetically, 
the change from subtraction to addition of pair separations in 
the BoP filter is equivalent to swapping an ((;, j), (k, I)) pair of 
ordered pairs to ((/, k), (j, I)). The filters are: 



(1) type I pair filter: 
\d(i, j) - d(k, 01 < e, 

(2) lifetime filter: 



|f,- - t k \ < At and \tj - f/| < At 
\tj - f/| < Af and \tj - t k \ < At 
otherwise: 



(3) 



X((i,j,k,l)) = 0, or 
X((i,j,k,[)) = l, or 
filter fails, $ x (4) 



(3) BoP filter: 



and 

and 

and 
or 

1 and 
and 
and 



\[x(i) - x(j)] - [x(k) - x(l)]\ < £BoP 

1 1X0 - 2/0')] - iy( k ) - y(0] I < fBoP 

\[z(i)-z(j)]-[z(k)-z(l)]\<e BoP ) 

\[x(i) - X(j)] + [ X (k) - X(l)]\ < £ BoP 

I [2/(0 - 2/0')] + ly(k) - 2/(0] I < £boP 
\[z(i)-z(j)] + [z(k)-z(l)]\<e B o P ). 

(5) 



Substituting a difference test for the BoP filter above, i.e. 
\d(i, j) - d(k, l)\ < £boP and \d(i, k) - d(j, l)\ < 6boP for the un- 
crossed and crossed pairs of pairs, respectively, would allow 
non-planar quadrilaterals (fold a rectangular sheet of paper along 
its diagonal to see this), whereas only a parallelogram can repre- 
sent a T 3 topological quadruple. 

A list of quadruples (pairs of pairs, each associated with a 
crossed quadruple parameter x) that satisfy all three successive 
filters is obtained. An arbitrary galaxy i is a member of Si > 
quadruples. The frequency of s in a given simulated set of 
galaxies is f(s), i.e. 



f(s) := |{i : s t 



(6) 



If a topological lensing signal is present, there should be a high 
number of galaxies which are a member of many quadruples, 
f(s) should be high for high s. This is a critical step in the 
successive filter method, introduced in Sect. 3.2.5 of Fujii & 
Yoshii ( 1201 II ): for s » 1, the statistic f{s) should be signifi- 
cantly higher in a catalogue containing topological lensing than 
in a simply connected catalogue. 



2.5. Statistical significance 

For an arbitrary realisation, let us define the cumulative number 
of galaxies that are each members of many quadruples 



(7) 



for some quadruplet membership number s* > 0, which for- 
malises the word "many" in "many quadruples". Thus, for a 
fixed value s*, F is a random variable which we model numeri- 
cally. To estimate what observational strategy is required to have 
a low chance of a false negative inference from the data, i.e. to 
find fi, the expected probability of falsely excluding the T 3 hy- 
pothesis, suppose that the observational result gives F - i for 
some non-negative integer i. Then the probability of falsely ex- 
cluding the Aurich (120081 ) T 3 hypothesis is the cumulative prob- 
ability P(F < i\T 3 ), which can be estimated numerically by find- 
ing the fraction of T 3 simulations for which F < i. Since we 
don't yet know the results of the observations, we have to weight 
this ("marginalise it") over the probability density function of i 
for the simply connected case. 

Thus, the expectation value of the probability of falsely re- 
jecting the hypothesis is 



Po 



£ (p(F = /|R 3 ) P(F<i\T 3 )). 



(8) 



i=0 



Similarly, to estimate what observational strategy is required to 
have a low chance of a false positive inference from the data, the 
expected probability of falsely supporting the T 3 hypothesis can 
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be written 

CO 

Y,(p(F = i\T 3 ) P(F>i\R 3 )). 



- s 



P(F = i\T 5 ) 



CO CO 



Y J p(F = j\R 3 ) 



2 2 = HT 3 ) {p(F 



1=0 /=! 

oo j 



J] J] [p(F = i\T 3 ) (p(F 

j=0 i=0 



J »] 



= Po, 



(9) 



i.e. the two expectation values are equal. 

For a given simulation parameter set, a set of simulations is 
analysed for each s* e {1, . . ., 10), a range that is found below 
(Sect. [3]) by inspection of N(s) versus s histograms for typical 
realisations. The minimum 



pi := minpo 



(10) 



determines the s* value for the optimal statistical test for the sim- 
ulation parameter set, and the expectation value of the probabil- 
ity of falsely excluding the T 3 hypothesis or of falsely detecting 
evidence in favour of it. 



3. Results 

3.1. Existing observations 

Table [5] lists observed objects and the implied positions where 
topologically lensed images should be observed if the T 3 hy- 
pothesis is correct, with zero uncertainty in the T 3 axis parame- 
ters given in TableQ}. Figure [3]shows the application of the suc- 
cessive filters method under the assumption that a survey of one 
square degree in the redshift range 5 < z < 6.2 is 100% complete 
in comparison to the surveys that found the high galactic latitude 
sample. The southern survey field is centred on the expected 
median celestial position of the implied objects (the fourth and 
fifth columns of Table [5). The topological signal is clearly very 
strong in Fig. [3] both for low Af = 0.01 hr l Gyr (top), in which 
case it is unlikely that stellar evolution would cause the earlier 
image of a galaxy to dim too much to be seen, and for high 
At - 1 h~ l Gyr, in which case dimming could weaken the test. 
The northern part of each catalogue analysed consists of the 44 
known objects, there are 41 topologically implied objects within 
the southern field, and the simulated southern subset consists of 
24 and 30 objects for Af = 0.01 h~ x Gyr and 1 h~ l Gyr), respec- 
tively. 

The uncertainty in the T 3 coordinates and fundamental 
length also need to be taken into account. An error of 2° in the 
vector (in the flat comoving covering space) from the observer to 
an object at a matched disc centre implies an error of about twice 
the size in the observing angle of the implied image, because the 
diameter of the shell containing objects at the distance of the 
matched disc centres is twice the radius of the shell. Thus, as 
Fig.|4]shows, a typical realisation of a 1 deg 2 southern survey is 
very unlikely to detect any topological quadruples (satisfying the 
filter criteria). The total disappearance of the topological signal 
in this figure (s — for all objects, since there are no topolog- 
ical images lying within the southern field at all), making the 
N{s) distribution even weaker than that of the simply connected 
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s (number of quadruples containing a galaxy) 



<H 10 




20 40 60 80 100 
s (number of quadruples containing a galaxy) 

Fig. 3. Quadruple frequencies N(s) obtained by the successive fil- 
ter method (Sect. \2.4\ applied to sets of objects in the multiply con- 
nected (dark curve, red online) and simply connected (pale curve, green 
online) cases, using the known, northern galactic objects in the SDF 
(second, third and sixth columns of Table |5}- The southern object 
set is calculated using the Aurich d2008t T 3 model in the multiply 
connected case, and simulated with an input auto-correlation function 
r = 10.0/T 1 Mpc,y = 1.8, r mm = 1.0/r 1 Mpc,r max = 100.0/r 1 Mpc in 
the simply connected case. The southern field is "observed" over 1 deg 2 
at the expected position in both cases. Top: Successive filter parameters 
(Sect. US are e = 1.0/T 1 Mpc.Af = 0.01/T 1 Gyr, e BoP = 1.0/T 1 Mpc. 
Bottom: Same, except that Af = IhT 1 Gyr. 



data set, is an artefact of the method. We did not add any "back- 
ground" simulated data around the northern galactic, observed 
SDF field, so the implied southern objects constitute a "survey" 
that contains less objects (zero) than the simulated southern sur- 
vey (16 objects in Fig.|4|i. Thus, let us consider fully simulated 
surveys. 

3.2. Planned or possible observations 

Given the uncertainties in the T 3 solution discussed in Sect. 12.11 
and listed in TableQ] let us consider a pair of northern and south- 
ern surveys centred on the matched disc centres for the same T 3 
axis, each of 64 deg 2 . As explained above, for a given known 
object, a two degree T 3 axis uncertainty yields an approximately 
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20 40 60 80 100 
s (number of quadruples containing a galaxy) 

Fig. 4. Quadruple frequencies N(s), as for Fig. [3] but a realisation of 
Gaussian uncertainties in the T parameters is applied, destroying the 
topological signal, since none of the implied topological images fall in 
the southern field in this realisation. 



four degree uncertainty for an observer placed halfway between 
the opposing images and searching for a topological image at 
the expected position. However, since we simulate both fields, 
other pairs of objects in the matched discs can be found, so 2° 
should approximately correspond to a 68% chance (for one an- 
gle offset from a Gaussian distribution of width 2°) of finding 
matching pairs. Thus, 64 deg 2 should approximately correspond 
to a 2cr (where <x is one standard deviation) chance of finding 
topological pairs. 

The uncertainty in L gives a lcr uncertainty of about 
150 hT x Mpc, i.e. a 2<x uncertainty of 300 h Mpc. For a sur- 
vey going from the matched disc centre to f3 ~ 4° away, Ta- 
ble |2] (for £2 m = 0.28) gives a typical redshift of z ~ 5.38, i.e. 
5.706 h" 1 Gpc in radial comoving distance from the observer. 
Inverting this, 5706 + 300 /T 1 Mpc gives 4.60 <z< 6.33. We 
extend this a little to be conservative, and increase the filter cri- 
terion e to 2 h 1 Mpc. A typical realisation is shown in the top 
panel of Fig. [5] Clearly, a pair of surveys that are wide enough 
on the sky and in redshift coverage, finding 300 objects above an 
idealised, z-independent completeness limit, has a good chance 
of strongly rejecting or supporting the T 3 hypothesis. 

Could photometric redshift estimates be sufficient? The 
lower panel of Fig. [5] shows the quadruple signal for a reali- 
sation with Gaussian peculiar velocity uncertainties simulated 
with cr(j6 pec ) = 0.005, i.e. a redshift uncertainty of <x(Az) « 0.03. 
The signal is clearly much weaker, but (in this realisation) is 
still distinguishable from the N(s) curve of the simply connected 
simulation. 

Figure [6] presents confidence levels, i.e. p\, as defined in 
Eq. ( [Tol l in Sect. 12.51 from ensembles of realisations. Unsurpris- 
ingly, low 0-(Az p hot) and high N g . d i generally give the statistically 
most significant results. The dependence on these two parame- 
ters is not fully monotonic. This is understandable because for 
a fixed cr(Az p h t), higher N„- d \ not only increases the numbers of 
topological quadruples, it also increases the number density, so 
that the numbers of non-topological quadruples also increase, 
yielding some complexity in the dependence. There is some sta- 
tistical noise in the contours. The total computing time for Fig. [6] 
on a single, recent 4-core processor would be from a few weeks 
to a few months (AppendixlBli. The actual calculations were per- 
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Fig. 5. As for Fig. [3] for fully simulated, multiply connected 
(dark, red online) and simply connected (pale, green online) catalogues, 
over northern and southern fields each of w 64 deg 2 centred at an- 
tipodal matched disc centres, with 4.3 < z < 6.6, and successive fil- 
ter parameters e = 2. Oh' 1 Mpc, At = l.0h~ l Gyr, e Bo p = 1/T 1 Mpc, 
and (in the simply connected case) an input auto-correlation function 
r = 10.0/r 1 Mpc, y = 1.8, r mm = 1.0/T 1 Mpc, r max = 100.0/T 1 Mpc 
[Eq. (|2)]. A total of 300 simulated objects (north plus south) are present 
in each catalogue. Top: Spectroscopic redshifts. Bottom: Photometric 
redshifts with cr(j8 pcc ) = 0.005, i.e. cr(Az) « 0.03. 



formed using parallel computing resources on several different 
machines. 

Figure [6] also shows that a typical photometric redshift er- 
ror of <x( Azphot) ~ 0.01 would enable rejection at p\ < 0.05 (i.e. a 
2cr rejection according to intuition for a Gaussian distribution), 
provided that A^ ga i > 500 in northern and southern surveys each 
over a> « 100 deg 2 . This is only a moderate confidence level 
for rejecting the hypothesis, but since p\ also represents the ex- 
pectation value of falsely accepting the hypothesis, it would be 
sufficient to provide a strong motivation for studying candidate 
topologically lensed galaxies further. A plot roughly similar to 
the lower panel of Fig. would be obtained, enabling the partic- 
ular galaxies most likely to be topologically lensed to be identi- 
fied. The particular realisation in the lower panel of Fig. |5]shows 
about 10 galaxies that are each members of s > 1 quadruples. 
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Fig. 6. Expectation values pi of falsely rejecting (or falsely accepting) 
the T 3 hypothesis for optimal choices of s* [Eq. UOH. as a function of 
photometric redshift < x(z ph ot ) [estimated as (1 + (z\ + Z2)/2)cr(B pc( .), using 
o"(/8 P cc) <k 1 and Eq. JA.4K where zi < z < Zi is the redshift range], and 
the number of galaxies N gal per combined (north plus south) catalogue. 
The lowest probabilities p\ (highest significance) results are for low 
°"(z P hot) and high A? ga i, at top-left in each of the panels. Top: 4.3 < z < 
6.6, solid angle per survey direction (south or north) a> « 100 deg 2 , 
Af = 1 lr l Gyr. The other parameters are identical to those for Fig.|5] 
Bottom: Same, except that At = 0.1 hT 1 Gyr. 



This is a small enough number for spectroscopic followup to be 
relatively easy to obtain. 

Comparison of the top and bottom panels in Fig. [6] illus- 
trates the potential role of evolutionary effects. The lower panel 
limits type II pairs to those with At = 0.1 hr x Gyr instead of 
Af = 1 h Gyr (upper panel). To attain moderate confidence, 
i.e. p\ < 0.05, photometric redshift errors would need to be tight- 
ened to about 0"(Az p hot) ~ 0.005, for slightly higher numbers of 
galaxies. For initial starburst durations not much shorter than a 
typical galaxy dynamical time of ~ 1 Gyr, Af = 0.1 hr x Gyr 
should give a high probability that a galaxy is included in a sur- 
vey at both its topological images. 

Figure [7] shows that increasing the survey area and redshift 
depth still further, to a> ~ 196 deg 2 and 4 < z < 7, would al- 
low statistically similar results for weaker accuracy of the pho- 
tometric redshifts and lower numbers of galaxies. An accuracy 
of cr(z phot ) < 0.02 or cr(z phot ) < 0.01, for Af = 1 h' 1 Gyr (upper 
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Fig. 7. Expectation values p x of falsely rejecting (or falsely accepting) 
the T 3 hypothesis, as for Fig. [6] but for a thicker, wider survey, over 
4 < z < 7 and uj » 196 deg 2 . Top: At = 1 hr 1 Gyr, Bottom: At = 
0.1 r 1 Gyr, 

panel) or Af = 0.1 lr l Gyr (lower panel), respectively, would 
give p\ < 0.05. 

Figure [7] also shows that for a sufficiently high number 
of galaxies, N ga \ > 900, a photometric redshift accuracy of 
cr(Azphot) ~ 0.01 would be sufficient for p\ < 0.01, i.e. an expec- 
tation value of false rejection of about 1%. Most of the present 
interest in photometric redshift techniques is at lower redshifts, 
but high statistical accuracy in z ~ 6 photometric redshifts would 
clearly be useful for deep redshift topological lensing. 

4. Conclusion 

Over the next few decades, wide-angle surveys in the z ~ 6 
redshift range will inevitably be performed. The results above 
indicate that with appropriate targetting and choices of observa- 
tional parameter limits, the speed with which the T 3 candidate 
for the topology of the Universe can be rejected or detected can 
be optimised. For more specific observational strategies for par- 
ticular telescope/instrument combinations, more detailed anal- 
yses could potentially be carried out by using detailed galaxy 
formation models, such as the hybrid -body/semi-analytic sim- 
ulations (Roukema et al. [T9931 [T9971 [200Tb that have been ex- 
tensively developed to simulate detailed galaxy properties as a 
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function of space and time (e.g., Hatton et al. 120031 and refer- 
ences thereof), including a specific focus on LAEs (e.g., Garel 
et al. 120121 1. More sophisticated combinations of existing ob- 
servations, including the VVDS 1400+05 field within 16° of the 
same southern axis and the VVDS 0226-04 field within 20° of 
the corresponding northern axis (Table |4} along with simulated 
future observations, should also potentially yield several alterna- 
tive strategies for topological lensing detection or rejection. 

The total number of L„ galaxies at 4 < z < 7, where L* is 
the characteristic luminosity of a Schechter function (Schechter 
11976b . is estim ated to be 0, ~ 3 x 10 3 /i 3 IVhxT 3 (e.g. Table 5, 
Bouwens et al. 1201 j} . This gives about 3-8 million L» galaxies 
for the pairs of 100 deg 2 and 196 deg 2 fields simulated above, 
over 4.3 < z < 6.6 and 4 < z < 7, respectively. For the mini- 
mum Afg a i Z 500 or Afg a i > 400 needed for achieving 5% expected 
confidence levels in the two field sizes, respectively (Sect. 13. 21 , 
surveys complete to L > 8.8L* or L > 9.9L*, respectively, would 
find the required numbers of galaxies. Because the computing 
time for the above calculations scales roughly as TV 4 , for suf- 
ficiently high iVg a i (since the heaviest computation is checking 
the list of possible quadruples; see AppendixlBl. simulations for 
surveys complete to L > L», i.e. for N gSL i > 10 6 , are clearly not 
practical without further filters added early in the successive fil- 
ter algorithm. 

Nevertheless, the most difficult aspect of predicting con- 
staints on the model is the difference between precise cosmol- 
ogy and accurate cosmology. The dark energy parameter Qa is 
suspected to be an artefact of using the FLRW metric (a homo- 
geneous solution of the Einstein equation) rather than the physi- 
cal average metric of the actually inhomogeneous Universe (e.g. 
Buchert & Carfora 120031 B uchert 120081 Celerier et al. 12010] 
Wiegand & Buchert 1201 Oi l. A significant confirmation of the 
T 3 hypothesis would provide a constraint on the inhomogeneous 
approach to cosmology, but this is unlikely to occur using an 
FLRW interpretation of the observational data if the time depen- 
dence of the FLRW metric parameters is too far from relativis- 
tically consistent formulae in the relevant redshift range (Larena 
et al. 120091 ). Some of the familiar FLRW relations may be al- 
gebraically valid in fully relativistic models, with an effective 
rather than a local physical interpretation (Buchert & Rasanen 
12012b . so detection of topological lensing might still be possible 
under the assumption of an FLRW metric. 
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Appendix A: Combination of cosmological and 
peculiar velocities in the FLRW model 



As derived in Synge (11960b and presented in Narlikar (11994) 
(see also Roukema 120 101 and refs therein), the expansion red- 
shift can be derived by parallel-transporting the four-velocity of 
the world line of a distant galaxy along a null geodesic (path of 
a photon) joining it to the observer. A fundamental observer (at 
rest with respect to the comoving spatial coordinate system) has 
zero peculiar velocity and a redshift of z cosm . The latter can be 
interpreted as a special-relativistic radial speed /3 cosm in natural 
units, where 
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Fig. A.l. Confidence level calculation time in seconds for quadruple 
searches of (the equivalent of) 30 realisations for a single parameter set, 
as a function of the total number of galaxies 2V ga i, for the simulations 
shown in Figs [6] and [7] for three different processor sets, as labelled. 
The times for the first processor set are halved to give the equivalent of 
four cores. Power-law fits to N gti \ > 400 are shown. 



where (f> cosm is the rapidity defined by J3 cosm 
ycosm := (1 -0t m T m , and rpec := (1 -A 2 



pec/ 



=: tanh0 cosm , 
y 1 * 2 . For a ra- 
dial peculiar velocity of y8 pec in natural units, the overall redshift 
Z can be calculated using Minkowski spacetime addition of the 
rapidities <p cosm and <p pec , where y8 pec =: tanh0 pec , i.e. the overall 
rapidity is 

4> = ^cosm + 0pec, (A.2) 

since addition of four-velocity vectors at the same spacetime lo- 
cation is meaningful. Thus, similarly to Eq. dA.lt . the overall 
redshift is given by 

1+z = e 

— g0cosm+0pec 

= (1 +Zcosm)(l +ZpecX (A.3) 



i.e. 



+ Z 



(A.4) 



Zpecl) ^ 1, this reduces to z — z 



When max(z c , 
first order in both redshifts. 

For high-redshift astronomy, i.e. z COS m > 1, the third term in 
the right-hand side of Eq. (IA.4t is more significant than the sec- 
ond, contrary to popular belief that sets the third term to zero. 
Nevertheless, for |z pec | ■« 1, which is realistic even for low- 
redshift astronomy, Eq. (IA.4t can be approximated to first order 
inAec, 



~l~ /^pec "t" Zqq&toP] 



pec ■ 



(A.5) 



We ignore gravitational redshift here, since observations close 
to the Schwarzschild radius are unlikely in the case of interest. 
Parallel transport of the four- velocity (see Narlikar [1994t implies 
analogous relations to those above. 

Appendix B: Simulation benchmarking 



(A.l) 



Figure IATI shows that for sufficiently high N ga \, the successive 
filter quadruple searches scale roughly as the fourth power of the 
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number of simulated galaxies. From top to bottom as labelled, 
the power law best fits are proportional to A^ 4 ?, N 4A ,, and A^ 3 ?, 

1 1 1 gal ' gal gal ' 

respectively. Parallelisation is via openmp. 



